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Abstract 



We present a kinetic theory for one-dimensional inhomogeneous systems with weak long-range 
1^ I interactions. Starting from the Klimontovich equation and using a quasilinear theory valid at order 

1/A^ in a proper thermodynamic limit N — > -|-oo, we obtain a closed kinetic equation describing 
the relaxation of the distribution function of the system as a whole due to resonances between 
different orbits. This equation is written in angle-action variables. It conserves mass and energy 
and increases the Boltzmann entropy (H-theorem). Using a thermal bath approximation, we derive 
a Fokker-Planck equation describing the relaxation of a test particle towards the Boltzmann dis- 
tribution under the combined effects of diffusion and friction. We mention some analogies with the 



'^ ' kinetic theory of point vortices in two-dimensional hydrodynamics. We also stress the limitations 

Q . of our approach and the connection with recent works. 
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I. INTRODUCTION 

Systems with long-range interactions have been the object of considerable interest in 
recent years [1]. Their dynamics is very rich and presents many interesting features [2]. 
Therefore, the construction of a kinetic theory adapted to such systems is valuable. Several 
kinetic theories have been developed in the past. The first kinetic theory was constructed 
by Boltzmann [3] in his theory of gases. In that case, the particles do not interact except 
during strong collisions. The extension of this kinetic theory to the case of particles in 
Coulombian or Newtonian interaction was considered by Landau [4] in the case of plasmas 
and by Chandrasekhar [5] in the case of stellar systems. Developements and improvements 
of the kinetic theory of Coulombian plasmas were made by Lenard [6] and Balescu [7] 
using more formal approaches allowing to take into account collective effects. They showed 
in particular how collective effects can regularize the logarithmic divergence of the diffusion 
coefficient at the Debye length. On the other hand, a kinetic theory was developed by Dubin 
& O'Neil [8, 9] and Chavanis [10, 11] to describe the "collisional" evolution of charged rods in 
a magnetized plasma or point vortices in two-dimensional hydrodynamics. In that case, the 
particles interact via a logarithmic potential in two dimensions. A specificity of this model 
is that the particles have no inertia so that the collision term acts in position space. More 
recently, Bouchet, Dauxois and Chavanis [12-16] have considered the kinetic theory of the 
Hamiltonian Mean Field (HMF) model, a toy model of systems with long-range interactions 
which possesses a lot of interesting properties and which can be studied in great detail. 
Finally, in very recent works, Chavanis [2, 17] and Benedetti et al. [18] have developed 
the kinetic theory of a 2D Coulombian plasma and Valageas [19] has developed the kinetic 
theory of a ID gravitational system in a cosmological context. 

Kinetic theories are important to describe the relaxation of a Hamiltonian system of par- 
ticles in interaction towards statistical equilibrium. They are also necessary to determine 
the timescale of the collisional relaxation and to prove whether the system will truly relax 
towards statistical equilibrium. Indeed, the convergence towards statistical equilibrium is 
not firmly established for complex systems with long-range interactions. A first reason is 
that systems with long-range interactions may exhibit non-markovian effects and spatial de- 
localization so that the Boltzmann if-theorem is difficult to prove and could even be wrong 
in a strict sense. Such problems have been discussed for example by Kandrup [20-22] in the 
case of self-gravitating systems and by Chavanis for point vortices [11] and for the HMF 
model [16]. These authors derive exact kinetic equations for these systems which apparently 
do not satisfy an if -theorem. The if -theorem is recovered only when additional assumptions 
(markovian approximation, local approximation, neglect of three body correlations,...) are 
implemented [34]. On the other hand, the relaxation of the system is due to a condition of 
resonance between different orbits and it may happen that the relaxation stops because there 
is no resonance anymore. This is the case for point vortices in two-dimensional hydrody- 
namics when the profile of angular velocity becomes monotonic [8, 11]. This is also the case 
for spatially homogeneous ID systems for which the Landau-Lenard-Balescu collision term 
cancels out [2, 13, 14]. In these cases, the system remains frozen in a long-lived "metastable" 
state (a stationary solution of the Vlasov equation) which is not the statistical equilibrium 
state. This means either that the system will not reach statistical equilibrium or that the 
kinetic theory is incomplete. Indeed, the collision term calculated in the above-mentioned 
kinetic theories corresponds to a term of order 1/A^ in a systematic expansion carried out 
in a proper thermodynamic limit A^ -^ +oo (one can rescale the parameters so that the 



coupling constant scales like k ~ 1/A^ while i? ~ A^, T ~ 1 and K ~ 1 [2]). Therefore, they 
describe the dynamics of the system on a timescale of order Nto where to is the dynamical 
time. If the collision term cancels out because of the absence of resonances, the relaxation 
towards statistical equilibrium (if any) will take place on a timescale larger than Ntn so that 
next order terms in the 1/A^ expansion must be taken into account. The calculation of these 
higher order terms has not been realized for any of the kinetic theories mentioned previously. 
However, for the HMF model, it is shown numerically that the relaxation towards statistical 
equilibrium takes place on a timescale of order N^'^to [23]. 

In all these works (except for point vortices), the system is assumed to be spatially 
homogeneous [17]. For the HMF model, a homogeneous phase exists if the energy is larger 
than a critical energy E^. For Coulombian plasmas, the interaction is shielded on a distance 
corresponding to the Debye length so that the system is spatially homogeneous. Stellar 
systems, on the other hand, are clearly inhomogeneous but most kinetic theories assume 
that the encounters can be treated with a local approximation [24] as if the system were 
spatially homogeneous. This results in a logarithmic divergence of the diffusion coefficient 
which is cured by introducing an appropriate cut-off at the Jeans scale. It is expected 
that this divergence would be removed if effects of spatial inhomogeneity were taken into 
account. Therefore, the description of the dynamics of inhomogeneous systems with long- 
range interactions, like the clustered phase of the HMF model below the critical energy Ec-, 
remains a challenging problem. In a previous paper [16], we have given some elements of 
kinetic theory for the inhomogeneous HMF model in a very special situation of low energy 
(close to the ground state) where the particles are strongly concentrated around a cluster 
at ^ = so that they all perform a harmonic motion with the same frequency. In that 
"pathological" case, some curious behaviors have been predicted. However, in more general 
situations, the particles will have different frequencies and we may expect that the relaxation 
will be due to some sort of resonances as in the case of 2D point vortices [8, 11]. 

In this paper, we present a kinetic theory for one-dimensional inhomogeneous systems 
with weak long-range interactions. Starting from the Klimontovich equation and using a 
quasilinear theory (Sec. II) valid at order 1/N in a proper thermodynamic limit N -^ +cxd, 
we obtain a closed kinetic equation describing the relaxation of the distribution function of 
the system as a whole due to resonances between different orbits (Sec. IV). This equation 
is written in angle-action variables (Sec. HI). It conserves mass and energy and increases 
the Boltzmann entropy so that a H-theorem can be proved (Sec. V). Using a thermal bath 
approximation (Sec. VII), we derive a Fokker-Planck equation that describes the relaxation 
of a test particle towards the Boltzmann distribution under the combined effects of diffusion 
and friction. Unfortunately, our derivation of the kinetic equation is based on a factorization 
hypothesis (28) that we cannot justify and that is probably not rigorously valid. However, 
we expect that more general approaches will be able to circumvent this difficulty. The 
important point is that the final kinetic equation (37) that we obtain is closed, without 
undetermined parameter, and that it possesses a lot of interesting properties (conservation 
laws and H-theorem). This gives some confidence in its relevance to describe the dynamics 
of the system even if the factorization hypothesis does not hold. In fact, in a very recent 
work, Valageas [19] has independently tackled a similar problem with a different method 
that avoids the factorization hypothesis. We shall discuss the connection between the two 
approaches in Sec. IX and show that our kinetic equation (37) can also be obtained from 
the approach of Valageas [19]. This provides therefore a rigorous justification for this kinetic 
equation without ad hoc assumption. 



II. THE KLIMONTOVICH EQUATION 

Let us consider a one-dimensional Hamiltonian system of particles in interaction described 
by the equations of motion 

dxi dH dvi dH 



dt dvi ' dt dxi ' 

1 ^ 
H = -^v1 + ^u{xi-Xj). (1) 

We assume that the particles interact via a weak long-range binary potential u = u{\xi — Xj\). 
This assumption is necessary to justify the quasilinear theory that we shall develop in the se- 
quel. This potential of interaction can be, for example, the potential u = —{k/27i) cos{6i—6j) 
of the HMF model [14]. We introduce the exact distribution function fd{x, v, t) = X]i=i S{x — 
Xi(t))S{v—Vi(t)) and the corresponding potential <l>(i(x,t) = f u{x—x')fd{x',v',t)dx'dv'. The 
exact distribution function satisfies the Klimontovich equation 

^ + [H„ M = 0, (2) 

where Ha = v'^/2 + ^a{x,t) is the exact energy of an individual particle and we have 
introduced the Poisson bracket 

dp dq dq dp'' 

with H = Hdi f = fd, q = X and p = v. The Klimontovich equation (2) contains exactly 
the same information as the Hamilton equations (1). 

We introduce a smooth distribution function / = (fd) which is the statistical average 
of the exact distribution function. We note $ = {^a) the corresponding potential. Then, 
we write fd = f + S f and $rf = $ + 5$. Inserting this decomposition in the Klimontovich 
equation (2), we obtain 

^ + ^ + [H,f] + [H, 5f] + [5$, /] + [5$, 5f] = 0, (4) 

where we have used H^ = H + 6H with H = v'^/2 + <l>(x, t) and 6H = 5$. Taking the 
average of Eq. (4), we get 

^ + [HJ] = -{[6<!>,6f]). (5) 

Subtracting this relation from Eq. (4) and neglecting nonlinear terms, we obtain an equation 
for the perturbation 

^ + [H,6f] + [5<^,f] = 0. (6) 

This linearization is the essence of the quasilinear theory [25] . This approximation can be 
shown to be equivalent to keeping only terms of order 1/A^ in a proper thermodynamic limit 
A^ -^ +00 [2]. In the A^ -^ +oo limit, the term in the right hand side of Eq. (5) vanishes 
and we obtain the Vlasov equation 

% + [H, f] = 0. (7) 



This equation admits an infinite number of stationary solutions of the form f = f{t) where 
e = f ^/2 + $(a:) is the individual energy. Such steady distributions can be reached through 
a process of violent relaxation involving only meanfield effects [26, 27]. If we keep terms of 
order 1/N or smaller, we shall obtain a kinetic equation taking into account the developement 
of correlations between particles due to granularities (finite N effects). This is similar to 
the effect of "collisions" in the theory of gases. Thus, the right hand side of Eq. (5) is 
similar to the collision term in the Boltzmann equation. This term is expected to drive the 
system towards statistical equilibrium. Therefore, it is expected to select the Boltzmann 
distribution among all possible steady solutions of the Vlasov equation (see, however, the 
discussion of Sec. V A). Since this term is of order 1/N or smaller, the effect of "collisions" 
is a slow process that takes place on a timescale Ntn or larger. Therefore, there is a scale 
separation between the dynamical time Id which is the timescale at which the system reaches 
a steady state of the Vlasov equation through violent relaxation due to inertial effects and the 
"coUisional" time tcou at which the system is expected to relax towards statistical equilibrium 
due to thermodynamical effects. Because of this scale separation, the distribution function 
will pass by a series of stationary solutions of the Vlasov equation depending only on energy, 
slowly evolving with time due to collisions (finite N effects). Hence, the distribution function 
can be approached by 

f{x,v,t)^f{e,t). (8) 

This means that the system is approximately in mechanical equilibrium [35] at each stage 
of the dynamics (due to violent relaxation) and that the "collisions" will drive it to thermo- 
dynamical equilibrium. The purpose of this paper is to obtain a closed expression for the 
collision term in energy space when the rapid dynamics has been averaged over the orbits. 
This is similar to the orbit-averaged- Fokker-Planck equation in astrophysics [24]. 

III. ANGLE- ACTION VARIABLES 

The previous decomposition has separated the smooth meanfield distribution /(e, t) which 
evolves slowly with time from the fluctuations 5f which are expected to have a decorrelation 
time on a much shorter timescale [36]. For inhomogeneous systems, it proves convenient 
to introduce angle-action variables [24] constructed from the slowly evolving distribution. 
Then, when we focus on the evolution of the perturbation, we can consider that the smooth 
distribution is "frozen". This will allow us to determine the collision term. Then, we can 
take into account the time dependence of the smooth distribution function through the 
kinetic equation (5). We define the action 

J = — ® vdx, (9) 

2-n J 

where 



v{x) = ^2{e-^x)). (10) 

By definition, the action is a function J = J(e) of the energy so it is constant on the 
iso-energetic lines in phase space. As discussed previously, we assume in a first step that 
the potential ^{x,t) is "frozen" on the timescale on which the fiuctuations have essential 
correlations so that J depends only on the energy (it is a new coordinate). We also introduce 
an angle variable 

f" dx 
<t> = n ^. (11) 



This angle parameterizes the iso-energetic hne with energy e. Therefore, a particle with 
coordinates (x, v) in phase-space can be described equivalently by the angle-action variables 
{(f), J). Noting that the period of an orbit with energy e is 

we find that Q = Q{e) = Q{J) is the pulsation of the orbit with energy e. We also note for 
future reference that J(e) represents the surface of phase space with energy smaller than e. 
Therefore, the surface of phase space with energy between e and e + de is g{t)dt with 

«W - f . (13) 

This is similar to the density of states in statistical mechanics (but here in the /i-space). 
Now, from Eq. (9), we have 



dJ 1 f dx 1 f dx 



de 271 J ^2(e - $) 271 J v Q{e)' 
so that g{e) = 1/Q{e). Accordingly, the density of particles with energy e is 



(14) 



m^) = mgie) = §^y (15) 

Since only the meanfield Hamiltonian appears in the advective term of Eqs. (5)-(6), it 
only involves the unperturbed motion of the particles. The smooth Hamiltonian equations 
for the conjugate variables {q,p) can be written 

dq dH dp dH 

dt dp ' dt dq' 

If we use the variables (a;, f ), we find that the dynamics is relatively complicated because the 
potential explicitly appears in the second equation. Therefore, this equation dv/dt = — $'(x) 
cannot be easily integrated except if $ = 0, i.e. for a homogeneous system. In that case, the 
unperturbed equations of motion reduce to x = vt + Xq, i.e. a motion at constant velocity. 
Now, the angle-action variables are constructed so that the Hamiltonian does not depend 
on the angles 0. Therefore, the Hamilton equations for the conjugate variables (0, J) are 

d(j) dH ^.^. dJ dH ^ ,^^, 

From these equations, we find that J = Cst. and = Q{J)t + 0o so that the equations 
of motion are very simple in such variables. They extend naturally the trajectories with 
constant velocity for homogeneous systems. This is why this choice of variables is relevant 
to develop the kinetic theory. Now, in terms of the angle-action variables, using the relations 

dHdSf dHdSf d6f 

^ '^^ 9 J 90 90 9 J d<p dJ' ^ ^ 



we find that Eqs. (5)-(6) take the form 

and 

IV. THE QUASILINEAR THEORY 

The equations (20)-(21) form the starting point of the quasihnear theory. The objective 
is to solve Eq. (21) for the fluctuation 6f and substitute the result back into Eq. (20). Here, 
we follow the approach of [30, 31]. We introduce the Fourier transforms [37]: 

(5/(0, J,t) = J2[ SfrrU uj)e'"''^e-'^'duj, (22) 

rn 

(5$(0, J,t) = J2 [ S^rnU uj)e'""^e-''''du;. (23) 

m 

Substituting these expressions in Eq. (21), we obtain 

df -^ 
[u - mn{J)] 6fm = -m-^5^„,. (24) 

Therefore, the perturbed distribution function is given by 

6f^{J,^) = ^-J^5^^{J,uj) + g^{J,uj), (25) 

u — m\l{J) 

where gm{J,uj) is solution to the equation 

[u-mn{J)]graiJ,io)=0, (26) 

and is related to the initial conditions [25, 30, 31]. Here, gm{J,uj) arises because of the 
discrete nature of the system, i.e. the fact that the exact distribution fd is a sum of 6- 
f unctions. 

The Fourier transforms of the potential and of the distribution function are related to 
each other by an equation of the form (Appendix A) 

5^^{J,UJ) =J2jAmm'{J,J')SU'{J',^)dJ', (27) 

m' 

where Amm'{J, J') are the Fourier components of the potential of interaction written in 
angle-action variables. At that stage, we shall make a rough approximation and assume 
that the function Amm'{J, J') can be factorized [38] according to 

A-mm'\J-, J ) = A.jn{J)A.m'{J )■ (28) 



Although this approximation is clearly a weakness of our approach, we can give the following 
arguments to motivate our study: 1. Without a decomposition of this type, it appears 
difficult to go any further in our approach. 2. There may exist some potentials of interaction 
for which this approximation holds exactly or approximately. 3. Such factorization could 
be obtained rigorously by using another system of coordinates or another basis on which 
the potential of interaction could be developed. 4. This approximation does not seem to be 
crucial: in particular the final kinetic equation (37) can be expressed in terms of Amm'{J, J') 
alone provided that we neglect collective effects. Therefore, the approximation (28) may 
lead to an inaccurate description of collective effects while preserving the correct structure 
of the kinetic equation. Assuming that Eq. (28) holds, we obtain 

m' 

This ratio is independent on m and J. Substituting Eq. (25) in Eq. (29) we find 

Am{J) 

Using the remark following Eq. (29), this can be rewritten 

1 + > ^ / Am'iJT "'^^f^,^ dJ'\ = V / A^iJ')grn'{J',uj)dJ'. (31) 



(30) 



AM) \'^^J ''"''^'> u;-m'Q{jT f ^ 



We introduce the dielectric function 

<^) = 1 + E / ^-'(^')^:7r5?T77A^^'' (32) 

so that 






S^UJ,^) = MZ) V j Arr,{J')grn'{J',Uj)dj'. (33) 

At this stage, we must require that e{uj) ^ for any real u. This is the case if the distribu- 
tion function /(e) is stable with respect to the Vlasov equation (see Appendix B). In fact, 
the distribution function /(e) which depends only on the energy results from a dynamical 
relaxation process as described at the end of Sec. II. This dynamical process (violent relax- 
ation) necessarily puts the system on a stable stationary solution of the Vlasov equation. 
Therefore, /(e, t) is a stable steady solution of the Vlasov equation at each stage of the 
dynamics which slowly evolves with time due to finite N effects (collisions). 
Substituting Eq. (33) in Eq. (25), we obtain 

SfmU^) = ""^^of n ^7^f ^ E j AM')9m'{J\oo)dJ' + gM,^)- (34) 

uj — m\l(J) e[ix)) ^-^ J 

m 

We can now compute the diffusion current using 



m.m' 



and the expressions (33) and (34). From that point, the calculations are similar to those 
carried out for the usual quasilinear theory [25, 30, 31] and they will not be repeated (see, 
e.g.. Appendix B of [2]). We finally obtain the kinetic equation 



dl 
dt 



ld_ 
2dJ 



m,m' 



(m'f](J'))|2 



mA„(J)M„/(/) 5(mfi(J) -m'fi(J')) 



x{f{J')m^-f{JW^Uf. 



(36) 



Using Eq. (28), the product Am{J)Am'{J') can be replaced by Amm'{J, J')- Then, we note 
that the function Am{J) occurs only in the dielectric function. From now on, we shall neglect 
collective effects and take \e(m'VL{J'))\ = 1. It is likely that our approach does not describe 
collective effects accurately. Therefore, the kinetic equation that we shall consider is 



Of _1 d y. 

'dt ~ 2dJ ^. 



mArnm'U Jr6{mn{J) - m'Q{J')) { f{J')m-^ - f{J)m' 



df 



dJ' 



df. 



(37) 
This equation depends only on Amm'{J, J') which can be calculated from Eq. (A5) for any 
potential. It does not depend on the individual terms Am{J) and Am' {J') which would be 
undetermined if Eq. (28) were not valid. Therefore, the kinetic equation (37) is well-posed 
for any potential. This may be an indication that its form does not crucially depend on the 
approximations that we have made in the course of the derivation. Furthermore, we shall 
see in the next section that this kinetic equation respects all the conservation laws of the 
Hamiltonian system (conservation of mass and energy) and increases the Boltzmann entropy. 
This again suggests that this equation provides a good description of the out-of-equilibrium 
dynamics even if our manner to derive it is not entirely satisfactory. This equation is the 
main result of the paper. It can be viewed as the counterpart of the Landau-Lenard-Balescu 
equation of plasma physics for spatially inhomogenous one-dimensional systems with weak 
long-range interactions. We note, in particular, that the relaxation is due to a condition of 
resonance mVL{J) = m'Q{J') between different orbits. This is similar to the kinetic equation 
obtained for point vortices in two-dimensional hydrodynamics [11]: 



dP 
'dt 



4r dr 



ridriS{Q{r) — ^{ri)) In 



1- — 



r dr ri dri 



(3J 



where Q{r) is the angular velocity, 7 is the circulation of a point vortex and r< (resp. r>) 
denotes the smallest (resp. largest) value of r and ri . Here again, the relaxation is due to a 
condition of resonance Q{r) = Q{ri) when the profile of angular velocity is non-monotonic. 
Equation (38) also neglects collective effects but the more accurate treatment of Dubin & 
O'Neil [8] allows to take them into account. 

V. PROPERTIES OF THE KINETIC EQUATION 



Let us derive some general properties of the kinetic equation (37). This equation can be 
written in the form of a conservative equation 



dt 



dQ 
' dJ' 



(39) 



with a current 



m,m' ^ '^ 

We impose the boundary conditions Q = at J = Jmin and J = Jmax so as to satisfy the 
conservation of mass. 



A. Equilibrium state 

The Boltzmann distribution 

/,, = Ae-^<'\ (41) 

is a stationary solution of the kinetic equation (37). Indeed, using Eq. (14), we have 



dfeq 



-Pfe,e\J) = -PfeMJ)- (42) 



dJ 
Introducing this relation in Eq. (40), we see that the integrand is proportional to 

{mn{J) - m'VL{J'))5{mVL{J) - m'fi(J')) = 0. (43) 

Therefore, the current vanishes and 

0. (44) 



dfeq 



dt 

An interesting question is whether the Boltzmann distribution is the only steady solution 
of the kinetic equation (37). Indeed, the current also vanishes if there is no resonance, i.e. 
if mQ{J) 7^ m'VL{J') for all orbits (m, J) 7^ (m', J'). In the case of point vortices described 
by a kinetic equation of the form (38), the Boltzmann distribution is a steady state but 
the relaxation stops as soon as the profile of angular velocity becomes monotonic, so that 
no resonance is possible, even if the system is not at statistical equilibrium. This is shown 
numerically in [33] . The question is whether a similar situation occurs in the more complex 
case of Eq. (37). Note, parenthetically, that we considered in [16] a case where there is no 
resonance. In that case, other terms must be taken into account in the kinetic theory and 
it is found that the system has a curious behavior. 

B. Conservation of mass and energy 

The conservation of mass 



M = jf{J,t)dJ, (45) 

is clear from the conservative form of Eq. (39) and the boundary conditions. Let us show 
that the kinetic equation (37) conserves the energy 

E= [ f{J,t)e{J)dJ. (46) 



10 



The time derivative of the energy can be written 

E= f^{J,t)e{J)dJ= fn{J)QdJ, (47) 

where we have integrated by parts and used e'(J) = Q{J). Then, introducing Eq. (40) in 
Eq. (47), we obtain 

E = -IY1 j dJdJ'QmA^m'iJ, J'f5{mn - rn'O.') [fm^ - fm'^\ . (48) 

m,m' ^ -^ 

Interchanging the dummy variables m, m' and J, J' we find that 

ra,m' 

(49) 

where we have used that Amm'{J, J') = Am'm.{J', J)- Summing these two expressions, we 
obtain 



E = — - y^ / dJdJ'Amm'iJ, J'Y{rnVL — mil')6{mil — mil') < f'f^-^ — f'^^'jrjj \ 



which is equal to zero according to Eq. (43). Therefore M = E = 0. 



(50) 



C. H-theorem 

Let us show that the entropy 

S=-j f{J,t)\nf{J,t)dJ, (51) 

increases. The time derivative of the entropy can be written 

Introducing Eq. (40) in Eq. (52), we obtain 

S=^ I dJdJ'-^ J2 mA^m'U J'f5{mn - m'n') Ifm^ " M'|^| • (53) 

•^ •' m,,m' ^ ^ 

Interchanging the dummy variables m, m' and J, J' we find that 

S = -^jdJdJ'^^ Yl ^'^^^'(J^ J'f6{mn - m'Q') h'm^ - fm'^\ . (54) 

m,m' ^ 

Summing these two expressions, we obtain 

S = - dJdJ'— ^ Ararn'{J, J'fS{mQ - m'Q') i fm— - fi^'-gj, ) ^ 



(55) 



which is positive. Therefore S* > 0. Note that the above properties persist for the kinetic 
equation (36). 

11 



VI. A MORE PRECISE KINETIC EQUATION 

In the preceding derivation, we have considered that the meanfield potential was "frozen" . 
In fact, this field evolves on a slow timescale and this evolution must be taken into account 
in the advective term of the kinetic equation (5). Let us write the kinetic equation in the 
form 

% + [H,f] = C{f), (56) 

where C{f) is the "collision" term. Since the distribution function (8) only depends on the 
energy (approximately), we have [H, f] ~ 0. On the other hand, 

^f(xvt)-^ + ^^ (57) 

where we have used de/dt = d^/dt since e = f ^/2 + $(a;, t). Therefore, the kinetic equation 
(56) becomes 

I- f !-(/)■ (-) 

The idea now is to average over the orbits in order to eliminate the dependence on (a:,f), 
using 

{X)ie,t) = ^-jf^—, (59) 

where we recall that §dx/v = 2Txg{e) is the density of states with energy e. The orbit- 
averaged-kinetic equation (58) then becomes 



1 f dx 
2Tr J ~ 

The first term is simply 



at + at ae ^^^^ 



1 rdxaf_aj af 
2^ J ~~dt ~ 'd^~at' 

Noting that 



0. (60) 

(61) 



a J 1/19$, 



the second term can be written 

1 rdxa^af_ aj af 

2^ J V~dt~de ~ ^~dt~de' 
Finally, the last term is equal to 



(63) 



1 /^C(/) = |^§ = ^. (64) 

2tt V (76 aJ ae 



Therefore, the final form of the kinetic equation is 



^^f{J')m^Jj-f{J)m'^}^dJ', (65) 



where, now, J = J{e,t). 

12 



VII. THE THERMAL BATH APPROXIMATION 

The kinetic equation (37) describes the evolution of the system "as a whole" . We have 
seen that it respects the microcanonical structure of the initial Hamiltonian system as it 
conserves mass and energy. Furthermore, due to the development of correlations between 
particles, the Boltzmann entropy increases. We shall now use this kinetic theory to describe 
the relaxation of a test particle in a thermal bath of field particles. Then, we interprete 
f{J,t) = P{J,t) in the kinetic equation (37) as the distribution function of the test particle 
and f{J',t) = fo{J') as the static distribution function of the field particles with which it 
interacts. This procedure transforms the integrodifferential equation (37) into a differential 
equation of the Fokker-Planck type [39] : 

^ = lijYl l^^rnn^'U J')H{^^{J) " m'f](j')) [W)^^ " ^('^)^'^} ^J' ■ 
m,m' 

(66) 
If we consider that the field particles are in thermal equilibrium, one has 

UJ') = Ae~^<''\ (67) 

Using 

^, = -pW)n{J'), (68) 

we obtain the Fokker-Planck equation 



dP 1 (9 /" r dP 1 

^ = 2 a7 5Z y ^^--'('^' J'f5{mn - m'Q')fo{J') Im— + pPm'Q'Uf. 

m,m' 



(69) 



Using the (5-function, this can be rewritten equivalently 



^ = lijJ2j ^--'('^' J'TKmn - m'n')m'fo{J') (^ + pPn) df. (70) 



Introducing the diffusion coefficient 



DiJ) = ^Y1 J^rnrn'iJ, J')H{mn - m'n')m' W)dJ\ (71) 



we obtain the Fokker-Planck equation 

dP _ d 
'dt~dJ 



D{J) (^^ + ppn{j)^ 



(72) 



This equation relaxes towards the Boltzmann distribution 

Pe{J) oc e'^'^^K (73) 

In terms of the energy, the Fokker-Planck equation (72) can be rewritten 






D{eMe) (^^ + PP 
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(74) 



Introducing N{e,t) = P{e,t)/Q{e), the density of particles with energy e, we obtain 



dN _ d 
~m~de 






where 



m,m' 

is the diffusion coefficient in energy space and 



(75) 



(76) 



V{e) 



6+ilnfi(e) 



is an effective potential. The stationary solution of Eq. (75) is 



N,(e) oc ^r^e"^^ oc g{e)e-'^% 



n{e) 



(77) 



(78) 



corresponding to the Boltzmann distribution. If we now consider a distribution of the field 
particles /o(^) such that mfi{J) ^ m'Q{J') for any (m, J) 7^ (m', J') so that this distribution 
is stable on a timescale Ntr, (see Sec. VA), the Fokker-Planck equation (66) becomes 



dP _ d 
'dt~dJ 



-wi|?--^ 



with a diffusion coefficient 



D{j) = bo{j)Yl 



ml 



WUT 



(79) 



(80) 



To arrive at this expression, we have used the fact that the only contribution to the collision 
term comes from m = m' and J = J' and we have simplified the (5-function according to 
6{m{n{J) - fi(J'))) = KJ - J')/\mn'{J)\. 

Note, for comparison, that if we consider the relaxation of a test vortex in a bath of field 
vortices with static vorticity distribution a;o(^), we obtain from Eq. (38) a Fokker-Planck 
equation of the form [10, 11, 33]: 



dP 

'dt 



7 d 
4r dr Jq 



r'dr'6{n{r) - n{r')) In 



r dr r' dr' 



If we assume that the field vortices are at statistical equilibrium ujq = N'-fPo = Ae ^'^^ 
(thermal bath), the Fokker-Planck equation becomes: 



dP _l d 

dt r dr 



\ or dr 



with a diffusion coefficient 



D(r) = - 



7 
4r2 



+00 



r dr'uJo{r')6{n - fi') In 



(82) 



(83) 
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To arrive at these expressions, we have used the (5- function and the fact that Q{r) = 
— {\/r)dil)/dr. Note that these expressions are vahd even if the profile of angular veloc- 
ity is non-monotonic since the Boltzmann distribution is always a stationary solution of the 
kinetic equation (38), and it is expected not to evolve at all (statistical equilibrium state). 
If we now consider a distribution of the field vortices ujQ{r) associated with a monotonic 
profile of angular velocity f2(r), so that it does not evolve under the effect of collisions on a 
timescale NId, we find that the relaxation of a test vortex is governed by a Fokker-Planck 
equation of the form: 



dP _l d 
dt r Or 

with diffusion coefficient 



rDir) I — P — ; 

or dr 



(84) 



D{r) = lj^^^\nNuoir), (85) 

where S(r) = rQ'{r) is the local shear created by the field particles (the angular velocity 
is related to the vorticity by ujq = (l/r)(r^f2)'). To arrive at this expression, we have used 
the fact that the only contribution to the collision term comes from r = r' and we have 
simplified the (5-function according to 6{Q{r) — Q{r')) = 6{r — r')/\Q'{r)\. The analogy with 
the previous equations can be noted. 



VIII. THE HOMOGENEOUS CASE 

In the homogeneous case where $ = and e = v^/2, we have 

(f) = x, J = v, n{J) = V. (86) 

On the other hand, 

where Um is the Fourier transform of the potential. Therefore, the kinetic equation (37) 
becomes 

dt 



2.^ A J2 I Hul6iv - v') (/'I - /|^) dv' = 0. (88) 

Thus, for a one-dimensional homogeneous system, the collision term vanishes [2, 13, 14] 
(this result is well-known in plasma physics; see last paragraph of [30]). As a result, the 
relaxation time is larger than NId. This has been shown numerically for the homogeneous 
HMF model where it is found that the system relaxes towards statistical equilibrium on a 
timescale A^^'^t/? [23]. By contrast, for inhomogeneous systems, the collision term in Eq. 
(37) does not necessarily vanish at order 1/A^ because there are new resonances. Therefore, 
we may expect that the system will relax towards statistical equilibrium on a timescale of 
order Nt£, provided that there are sufficient resonances to drive the relaxation to completion 
(see the remarks in Sec. VA). 

On the other hand, for a one-dimensional homogeneous system, the velocity distribution 
of a test particle evolving in a bath of field particles with static distribution fo{v) is governed 
by the Fokker-Planck equation 
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If the field particles are at statistical equilibrium with the Maxwellian distribution function 



/o oc e" 



-fiv^l2 



(thermal bath), the Fokker-Planck equation becomes 



dP 

'dt 



d_ 
dv 



dP 

D(v){— + l3Pv 



with a diffusion coefficient 



D{v)=2n'fo{v)Y, 



ul^\m\ 



(90) 



(91) 



proportional to the distribution function of the thermal bath. In fact, since the kinetic 
equation vanishes at order 1/A^, we can consider a bath with any stable distribution of the 
Vlasov equation; it will not change on a timescale Ntu. On the other hand, collective effects 
can be included in the kinetic theory through the dielectric function. These generalizations 
are discussed in [2, 13, 14]. They lead to a Fokker-Planck equation of the form 



with a diffusion coefficient 



dp d 

dt dv 


\Div)(f-Py^)] 
\dv dv J 


D{v) = 


"2 1 

^^^ e(m,mv) 



(92) 



(93) 



where the dielectric function is 



e(m, cuj 



1 + 27r-u^ 



f» 



ijj jm — V 



dv. 



(94) 



These kinetic equations are again very similar to the previous equations. For a Maxwellian 
distribution of the field particles /o = {(3/2'kY^'^p exp(— /9f^/2) (thermal bath), the dielectric 
function can be explicited [2]. The Fokker-Planck equation is then given by Eq. (90) with 
a diffusion coefficient 



D{v) = 2rc'fo{v) Y, 



\m\ ut 



l + 27rUmPpB{Jllv 



(95) 



+ 27r3/33p2^2^y2e-/3«2 



where B{x) = 1 — 2xe ^ L e* dt. The result (91) is recovered if we neglect collective effects 
and take e{m,mv) = 1. 



IX. CONNECTION WITH OTHER WORKS 

While this paper was in course of redaction, we came accross the very interesting study 
of Valageas [19] for ID self-gravitating systems. This author develops a kinetic theory for 
inhomogeneous systems by starting directly from the Hamiltonian equations of motion and 
using a perturbative expansion of the trajectories in powers of 1/A^. He can therefore obtain 
the coefficients of diffusion and friction directly from the equations of motion. Although the 
methodology is different, the expressions of the diffusion and friction coefficients obtained 
by Valageas [19] are very similar to those obtained in our approach (we shall see, in fact. 
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that they coincide after simple transformations). Since this author does not make the 
factorization hypothesis (28) this shows that this assumption can be avoided by using a 
different approach. We note, however, that the approach of Valageas does not take into 
account collective processes encapsulated in the dielectric function (polarization cloud). This 
is precisely in order to obtain this function that the factorization (28) was necessary in our 
approach. Therefore, the agreement between the two approaches when collective effects are 
neglected is consistent. 

There are, on the other hand, differences between the two approaches because Valageas 
[19] finally obtains a Fokker-Planck equation of the form 



f = S*^-'^' + 1'^"^^*- (^«' 



leading to a stationary solution 



Neie) oc 7^7-e"'^^ (97) 

DJe) 



which differs from the expected Boltzmann distribution (78). In fact, in order to arrive at 
Eq. (96), Valageas makes some approximations that are, in fact, not necessary. Indeed, he 
obtains the following expressions for the diffusion and friction coefficients 






(9J 



/ (A J)^ 



. ^^ = I dJ'fo{J')Y,^rnm'{J,J'fm'5{m^-m'Q.'), (99) 

m,m' 

where we have used our notations to facilitate the comparison with the equations of the 
present paper. If we assume that the dynamics can be described by a Fokker-Planck equa- 
tion, the evolution of the distribution function is given by 



dP^l^ (/{AJl\ \_^ // AJ\ 
dt 2dJ^ W At / J dJ \\At / 



(100) 



The Fokker-Planck equation (100) can be rewritten 



dP d f dP \ , , 

where we have introduced the diffusion coefficient 

and the friction term 

Now, using Eqs. (98) and (99), we obtain 

D = ]^ fdJ'foiJ') J2 Aarn'U J'fmH{m^ - m'fi'), (104) 

m,m' 
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and 

V = -\j dJ'U{J')Y,^'^A^m'{J,J'fm5{mQ-m'Q'). (105) 

m,m' 

Integrating by parts the second equation, we have equivalently 

Tj = — I dJ'——y Amm'i-J, ■J'Y'm''m.'5{mi^ — m'i^'). (106) 



- / dJ'—— y^ Ammi{J, J')'^mm'5{mil — mil' 

m,m' 



Therefore, the Fokker-Planck equation (101) can be rewritten as 
dP _1 d ^^ fmA^^,{J, J'f5{mn{J) - m'n{J')) S fo(J')m^ - P(J)m'^ \ dJ', 



m,m' 



— — -^-^ /_^ I iii-^mm'y^-, o ) uyiiL^Lyo ) — nu ^uyo )) \ joyo J"''-KJ ~ ^ \o jnt — 



(107) 
which coincides with our kinetic equation Eq. (66). Therefore, the two theories give the same 
results and the expressions of the diffusion and friction terms (98)-(99) can be obtained from 
our approach by starting from Eq. (66) and using the steps (98)-(107) in the other way round 
as done in [2]. When we perform a thermal bath approximation as in Sec. VII, we obtain 
from Eq. (106) that 7] = —Dj3Q{J) which can be interpreted as an appropriate Einstein 
relation. Substituting this relation in Eq. (101), we obtain Eqs. (72) and (75) without 
approximation. These Fokker-Planck equations relax towards the Boltzmann distribution 
(78), contrary to Eq. (96) obtained by Valageas [19]. This is more satisfactory on a physical 
point of view. Then, if we account for the fact that the distribution of the bath in Eq. 
(107) is not stationary but evolves slowly under the effect of "collisions" in a self-consistent 
way, we obtain the integrodifferential kinetic equation (37) or, more precisely, Eq. (65). 
Therefore, we can start from the approach of Valageas (which does not make the factorization 
hypothesis (28)) to rigorously justify the kinetic equation (37). This equation was not given 
by Valageas [19]. It describes the evolution of the system as a whole and conserves mass and 
energy. As we have seen, these conservation laws and the H-theorem directly result from 
the symmetrical form of the kinetic equation (37). 

X. CONCLUSION 

In this paper, we have introduced a kinetic equation (37) or (65) that describes the 
dynamics of spatially inhomogeneous one-dimensional systems with weak long-range inter- 
actions in angle-action variables. We have shown that this kinetic equation conserves mass 
and energy and increases the Boltzmann entropy (H-theorem). The evolution of the system 
is due to resonances between different orbits. The Boltzmann distribution is a stationary 
solution of this equation, but it is not necessarily the only one. The system can be frozen in 
another steady state of the Vlasov equation if there is no resonance, i.e. if mVL{J) ^ m'Q[J') 
for all orbits (m, J) 7^ (m', J'). On the other hand, using a thermal bath approximation, we 
have obtained a Fokker-Planck equation (72) that describes the relaxation of a test particle 
towards statistical equilibrium. These nice properties suggest that these kinetic equations 
are relevant to describe the evolution of the system. Unfortunately, the derivation of the 
kinetic equation (37) relies on a factorization hypothesis (28) of the Fourier components of 
the potential in angle-action variables which is not rigorously justified. However, we have 
shown that the recent work of Valageas [19], which uses a different approach where this 
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assumption is not required, leads to the same results after some transformations. Therefore, 
we can start from this approach to rigorously derive Eqs. (37), (65) and (72) which were 
not given in [19]. Thus, these equations should provide a good description of the dynamics 
of inhomogeneous systems with long-range interactions. These results can be applied to the 
HMF model and to the dynamics of point vortices in two dimensions. This will be presented 
in a future work [33]. Other extensions are also possible. 
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APPENDIX A: THE FUNCTION A„„ ( J, J') 

The relation between the potential and the distribution function is given by 

$(a;) = f u{x-x')f{x',v')dx'dv'. (Al) 

Introducing angle-action variables such that x = x{(f), J) and v = v{(f), J) we obtain 

$(0, J)= f u(x - x')f{(p', J')d(p'dJ', (A2) 

where we have used dxdv = d(f)dJ. Taking the inverse Fourier transform of Eq. (23) and 
using Eq. (A2), we find that 

^M) = [u{x- x')f{<P', J')e-^-^'^^d<P'dJ'. (A3) 

Inserting Eq. (22) in Eq. (A3) we get 

^M) = Y.ju{x- x')LiJ')e^^"^'*'-"^'^^^d<P'dJ'. (A4) 

m' 

Defining 

A„.m'U J') = ^ /^(^ - x')e^(™'<^'-™^)rf0#', (A5) 



we finally obtain 



(A6) 



APPENDIX B: DISPERSION RELATION FOR THE VLASOV EQUATION 

Let us consider the linear dynamical stability of a stationary solution /q = fo{J) of the 
Vlasov equation 

% + [H, f] = 0. (Bl) 
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If we substitute the decomposition f = fo + S f in the Vlasov equation (Bl) and neglect 
nonhnear terms, we find that the perturbation 6f obeys the hnearized Vlasov equation 

The Vlasov and linearized Vlasov equations (B1)-(B2) are similar to the Klimontovich and 
linearized Klimontovich equations (2)-(21). It is important to note, however, that in the 
Vlasov equation the distribution functions / and 5f are smooth distributions while in the 
Klimontovich equation they are expressed in terms of (5-functions. Writing the perturbation 
in the form 

<5/(0, J, t) = ^ 5/„( J, cu)e^'"^e-^-*, (B3) 

m 

we can proceed exactly like in Sec. IV except that (7 = since the perturbation 5f of the 
smooth distribution / does not involve any (5-function. Therefore, we find that e{uo) = 0, i.e. 



e{uj) = 1 + > / A^iJ'Y '-^^kr^df = 0. (B4) 



^J LO-m'Q{J') 

This is the required dispersion relation from which we can study the linear dynamical sta- 
bility of the steady state /o = fo{J) with respect to the Vlasov equation [32]. Note that 
the integration must be done by using the Landau contour. In general, u = Ur + ioJi is a 
complex pulsation, except at the point of marginal stability cUj = where it is real. In the 
kinetic theory of Sec. IV, we must require that the smooth meanfield distribution /(e) is 
Vlasov stable, which is the case on the basis of physical grounds (see discussion in Sec. II). 
In that case, e{(jJr) 7^ for any real uir- 
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